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Abstract 

All spacetimes for an irrotational collisionless fluid with a purely 
electric Weyl tensor, with spacetime curvature determined by the ex- 
act Einstein field equations, are shown to be integrable. These solu- 
tions include the relativistic generalisations of the Zeldovich solutions 
of Newtonian theory. Thus our result shows the consistency of vari- 
ous studies of " silent" universes (where such consistency was assumed 
rather than proved). 



Subject headings: 
cosmology - galaxies: clustering, 
formation - hydrodynamics - relativity - exact solutions 



1 Introduction 



The dynamics of self- gravitating collisionless matter (dust) that is irrota- 
tional and with vanishing magnetic part of the Weyl tensor has been the 
focus of interest in the last couple of years || 0, 0. 

The assumption that the matter can be described by a collisionless fluid 
is reasonable since we wish to study cosmological perturbations during the 
matter dominated era. However, because of the high degree of non-linearity 
in Einstein's field equations, a number of additional physically motivated ap- 
proximations have been introduced to make the problem more tractable. 

The second assumption is that the fluid flow is irrotational. Physically 
this is reasonably well justified provided the scales of interest are not too small 
(i.e. where rotation starts to become important). Then the kinematical con- 
dition that uj a b = implies the flow is hypersurface orthogonal. Furthermore 
if it is initially irrotational then it will remain irrotational in the absence of 
dissipative effects. It is the third assumption of vanishing "magnetic" part 
of the Weyl tensor that has led to most discussion. 

A remarkable feature of perfect fluid spacetimes with vanishing vortic- 
ity ujab and "magnetic" part of the Weyl tensor H a b is that they admit an 
orthonormal tetrad associated with the matter 4- velocity u a which is si- 
multaneously an eigen- tetrad for the shear of the matter flow o a b and the 
"electric" part of the Weyl tensor E ab [|I|. It follows that, apart from two spe- 
cial cases, coordinates exist in which the metric g a b and the tensors a a b and 
E a b are all diagonal. This together with the further condition that the flow 
is geodesic ii a = leads to a considerable simplification of the propagation 
equations, reducing them to a set of six ordinary differential equations (the 
Raychaudhuri and continuity equations together with two pairs of equations 
for the independent components of a a b and E a b) ||. In addition to these 
propagation equations various constraint equations have to be satisfied that 
result from these conditions. Apart from the equations that directly express 
the restriction H a b = 0, there is an additional constraint which follows from 
the H Bianchi identity. Physically this can be thought of as limiting the 
spatial variation of the "electric" tidal field E a b. 
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When these constraint equations are satisfied, each fluid element evolves 
independently of each other until the formation of caustics, when the one- 
to-one mapping between fluid elements and space points is lost. Such a 
universe model was dubbed silent by Matarrese et. al. |7| , because no infor- 
mation can be exchanged between the fluid elements which was not already 
present in the initial conditions (which must be chosen to satisfy the con- 
straint equations). 

The aim of this paper is to show that the required integrability condi- 
tions are satisfied: the time-propagation equations are consistent with the 
constraints. Two of the constraints are clearly sensitive to the assumption 
of Hij = and warrant a full consistency check: viz i) the usual H = = w 
constraint which restricts the spatial distribution of the shear a^; ii) and the 
additional constraint due to = described above. In fact in showing 
that in the case of dust, vectors of the tetrad are hypersurface orthogonal, 
Barnes and Rowlingson (1989) 0] performed a consistency test on the part 
of the constraint for which i — j. We carry out a full consistency on all 
the constraints. 

Sections |3| and |] list the evolution and the constraints equations in co- 
variant form. Results of calculation of the time derivative of constraints are 
given in section |5|, using a tetrad formalism. A sample calculation appears 
in the appendix. Consistency analysis is performed first for the trivial Type 
O fields, which evolve as Friedmann-Robertson- Walker (FRW) models. The 
type D and type I fields are tested separately. In both cases, integrability 
conditions are consistently satisfied. 

Kinematic and tetrad notation used here are the same as in Ellis (1967) 
||. Latin indices run from to 3 and Greek indices from 1 to 3. The tetrad 
{eo, e^} is chosen with the timelike vector eo chosen as the fluid flow vector, 
and the associated parameter r measuring proper time between the surfaces 
orthogonal to the fluid flow (spanned by the vectors ej}). This can be done 
since for dust we have zero acceleration. Such a specialization simplifies time 
propagation calculations. 
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2 Relativistic dynamics of irrotational dust 



In this section we will give a brief summary of the relativistic dynamics 
of irrotational dust [3] in terms of covariant variables that represent the 
observable kinematical and dynamical quantities [4], focusing on the case of 
irrotational dust with vanishing "magnetic part" of the Weyl tensor: p = 

Uab = H ab = 0. 

2.1 Variables 

The relativistic dynamics of dust is determined by the Einstein Field equa- 
tions and by the continuity equation for the matter stress-energy-tensor 
Tab — pu a Ub, where p is the energy density and u a is the normalized 4- 
velocity of the fluid (u a u a = — 1). At each spacetime point we can define 
a projection tensor h ab = g a b + u a u b for which h ab u a = 0. With u a and h ab 
it is possible to split the covariant derivative of any tensorial quantity into 
a time derivative and a spatial derivative. In particular the first covariant 
derivative of u a can be written as u a -b = v ab + a a u b , where v ab = h c a h d b u c -d is 
the spatial part satisfying v ab u h = 0, and a a = u a = u a ;b u b is the acceleration 
vector which is orthogonal to the fluid flow (a a u a = 0). It is standard to 
split Vab m to three kinematical quantities: it's trace = v a . a , symmetric, 
trace -free part a ab = — ^h ab Q, and antisymmetric part u ab = V[ ab ]. The 
physical meaning of these quantities are as follows: 9 measures the volume 
expansion or contraction of the fluid, a ab is the shear tensor which describes 
the rate at which a spherical fluid volume is distorted into an ellipse and 
u a b is the vorticity tensor giving the rate of rotation with respect to a local 
inert ial frame. 

2.2 Kinematic equations 

For a perfect fluid the fluid acceleration is only determined by pressure gradi- 
ents so the restriction of vanishing pressure implies that a a = 0. This means 
that each fluid element moves along a geodesic. The conservation of energy 
and momentum T ab . b = leads to only one further equation, the continuity 
equation: 

p = - P Q. (1) 
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With the second restriction of vanishing vorticity io ab = 0, the equations for 
the kinematic quantities follow from the Ricci identity: u a -d- c —u a - c .d = R a bcdU b 
[4,5]. The expansion scalar 6 obeys the Raychaudhuri equation: 

6 + |0 2 + 2a 2 + \k P = , (2) 

where a 2 = ^cr ab a ab is the shear scalar and k = 8itG is the gravitational 
constant. Note that in the homogeneous and isotropic case where a ab = 
and 6 = 3H, this equation reduces to the familiar Friedmann equation: 

3H + 3H 2 + \k P = . (3) 

The only other kinematic equation is for the shear. It is given by: 

<j a b + 0~acO~ C b - l& 2 Kb + |0O"a6 + E ab = , (4) 

where E ac = Ei ac \ = C a b c dU b u d is the "electric" part of the Weyl tensor 
Cabcd (satisfying E ac u c = 0, E a a = 0). E ab is that part of the grav- 
itational field which describes tidal interactions. The Weyl tensor can 
be decomposed into E ab and another tensor called the "magnetic" part: 
H ac = H {ac) = \ri ah 9h C ghcd u b u d (satisfying H ac u c = 0, H a a = 0). This is 
the part of the gravitational field that describes magneto-gravitic effects, 
and allows gravitational waves. 

So far we have explicitly made two assumptions, those of vanishing pres- 
sure and vorticity. In what follows we will make the third and most important 
assumption, namely, neglecting the influence of the "magnetic" part of the 
Weyl tensor. This implies neglecting the interaction of gravitational waves 
with the system. 

2.3 Constraint equations 

Besides the evolution equations for the kinematical equations, there are sev- 
eral constraints that our variables must satisfy. On setting p = u ab = H ab = 
we obtain as non-trivial constraints, 

h e b (29* _ h d c(T bcj = o ? (5) 
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the '(0, z/)' field equations, and 



tfahUvfyjfbcU* = , (6) 

the condition that H a b = 0. Additionally there is a Friedmann like equation 
giving the Ricci-scalar of the 3-spaces orthogonal to u a ; this is a first-integral 
of (1) and (2). 

2.4 Bianchi Identities 

Additionally, the Bianchi identities must be satisfied, as they are the inte- 
grability conditions for the other equations. With our restrictions, they take 
the form: 

h\h d s E as , d = \h\p* , (7) 

h m a h t c E (M: + h mt a ab E ab + QE rnt - 3E/ n a t)s = -\pa tm , (8) 
the 'div E 1 and equations respectively, and 

rf^u b a d p E qd = , (9) 

and 

hJ m V t > sd u r E a s . 4 = . (10) 

which are the 'div if' and 'if' equations respectively in the case where H ab 
vanishes. 

3 Tetrad approach 

A direct conversion of these equations into a tetrad system which is an eigen- 
frame for both the shear tensor and the Weyl tensor yields the following 
time-evolution equations (Barnes and Rowlingson, 1989 [|l|]); 

P = -po , 

9 = -|0 2 - W + *2 2 + <*») - \ P , 

cr» = -K) 2 -|^ M + |(o r i 2 + (T2 2 + a3 2 )-^, 

4 = -OEp - Ipa^ + 3a^ - (a l E 1 + a 2 E 2 + a 3 E 3 ) , (11) 
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being respectively the tetrad forms of (1), (2), and of the diagonal parts of (4) 
and (8). The non-diagonal parts of (4) and (8) (the propagation equations 
for E ab and & a b with a ^ b) introduces two additional conditions; 

o = K-^rv. (12) 

= (Ev-EJVon fx^u, (13) 

where the Ricci rotation coefficients are defined by T a b c := e a . Vfc e c [3], and 
one can raise and lower indices in the first place as usual by use of the metric. 

3.1 Uniqueness of tetrad 

The tetrad used here is an orthonormal tetrad. If we denote the local coor- 
dinate system by {x 1 } and the tetrad by {e a } then the equations (Ellis, 1967 

I) 

e a = ei{d/dj) (14) 

define the functions e l a , which are components of the tetrad vectors e a with 
respect to the basis d/dx\ and are also directional derivatives of the coordi- 
nate functions x % as 

e\ = d a (x l ). (15) 

The choice of a tetrad which simultaneously diagonalizes the shear tensor 
and the Weyl tensor possess the following uniqueness properties: 

1. The tetrad is uniquely determined for the two cases: 

(a) Distinct and/or distinct values. 

(b) Degenerate shear tensor and Weyl tensor, with the principal 
planes not coinciding i.e., (say) E\ = E 2 ^ E 3 and o\ = cr 3 ^ a 2 . 
Here equations (|12"D and flT3|) yield 

PV = 0; fi ^ v . (16) 

2. The tetrad is free by a rotation in (say) the e\,ei plane if both the 
shear and the Weyl tensor are degenerate in the same plane i.e., E\ = 
E 2 ^ E 3 and 01 = a 2 ^ cr 3 . 
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We write the tetrad freedom as 

e\ = e\ cos + e 2 sin , 

e 2 = — e\ sin <f> + e 2 cos 4>'i , 4> = 4>(x°,x l ) (17) 
and from ( |T2]) or (|H|) we have 

r 3 02 = r 3 i = o . (is) 

A rotation of the tetrad can then be performed to determine d<p/dx° 
from the requirement that 

r 2 i = o . (19) 

This leaves the tetrad arbitrary by an initial rotation in a surface x° = 
const. The remaining freedom is later used to set r 1 32 = in ([51]) 
below. 

Equation (|i~6|) is thus valid generally and simplify calculations greatly. Both 
equations fll2[) and (|13|) are satisfied. Also it follows trivially from the prop- 
agation equations that E^ = if and only if <r M = (FRW model). 



4 Constraint Equations 

The constraint equations give spatial restrictions on the tetrad form of the 
dynamical variables. The (0, v) field equations (5) are: 



w = 


0101 + (0i 


-02 


)r 2 2 i + (01 


-0-3! 


r 3 3i , 


(20) 


w = 


<9 2 02 + (02 


-01 


)T\2 + (02 


-03; 


r 3 32 , 


(21) 


PsO = 


^30"3 + (0"3 


- 01 


(r 1 ^ + 


- 2 ; 


r 2 23 . 


(22) 



The 'div E 1 equations (7) take the form 



1 Ap 


= d 1 E 1 + (E 1 


— E 2 


)r 2 2i 


+ (Ex 


— E 3 


)r 3 3i , 


(23) 


\d 2P 


= d 2 E 2 + (E 2 


-Eh 


)T\2 


+ (E 2 


-E 3 


ir 3 32 , 


(24) 




= d 3 E 3 + (E 3 


-E x 


)T\ 3 


+ (E 3 


-E 2 


>r 2 2 3 • 


(25) 



The 'div if' equations (9) are identically satisfied by the tetrad variables 
(indeed it is they that allowed us to simultaneously diagonalise the shear and 
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the Electric part of the Weyl tensor). 

The equation (6) for vanishing takes the form 



V x 32 {a 2 - a x ) = r\ 3 (a 3 -at), (26) 

r 2 3 i(a!-a 2 ) = r 2 13 (a 3 - <x 2 ) , (27) 

T\ 2 (a 2 - a 3 ) = T^iax-as), (28) 

d x {a 3 -a 2 ) = r 3 3 i(a 1 -a 3 )-T 2 21 (a 1 -a 2 ) , (29) 

d 2 {a 3 -a 1 ) = r 3 3 2 (o-2 - o- 3 ) - r 1 i 2 (cr 2 - o"i) , (30) 

d 3 (a 2 -a 1 ) = r 2 23 (0-3 - cr 2 ) - r\ 3 (0-3 - a t ) . (31) 

Finally the dynamical restriction = 0, equation (10), introduces con- 
straints given by 



(32) 
(33) 
(34) 

(35) 
(36) 
(37) 



5 Time propagation of constraints 

For the propagation equations to be integrable with the chosen restrictions on 





r 32 (e 2 — 


Et) = 


r 1 2 3(-E3 — Ei) , 






r 2 3i(-E'i — 


E 2 ) = 


r 2 i3(-^3 — e 2 ) , 






r 3 i 2 (£ , 2 — 


E 3 ) = 


r 3 2 i(-Ei — e 3 ) , 




di(E 3 


-e 2 ) = 


r 3 3i(-£'i 


— E 3 ) — T 2 2i(Ei 


-E 2 


d 2 (E 3 


-£1) = 


r 3 32 (s 2 


— E 3 ) — T 1 12 (£' 2 


-Ei 


d 3 (E 2 


-#1) = 


r 2 2 3(£ , 3 


— e 2 ) — y i 3 {e 3 


-Ei 



the kinematic variables (w^ = = H^), the constraint equations (|2Q|- p 7 ^) 
must be preserved during the time development of the system. We focus here 
on the H constraint (32-37); a similar analysis is also valid for constraints 
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The time propagation of the H constraint in covariant form (10) gives 

= hKrf )klm {(E n i);m - E\ p vP, m - R n mpq ESu p + W mpl E n q u p ] . (38) 

If we convert the new constraint (|38| ) to the tetrad form we obtain the fol- 
lowing. For /i = v, we find 

(39) 
(40) 
(41) 

which can be written equivalently as: 

= r 1 32 (E 2 -E 1 )(a 2 -(7 3 ) , (42) 

= r 2 13 (^3- ^2)^1-^3) , (43) 

= r 3 21 (£i-£3)(<7i-<7 2 ) ■ (44) 

For fi 7^ v, we get 






= r 1 23 (-E 3 


-E l )(a 2 


- ^3 





= r 2 3 i(£i 


~E 2 )(a 1 


- ^3 





= r 3 12 (-E 2 


-£ 3 )(*i 


- <T 2 



= (E 2 -E 3 )d 1 a 3 + (a 2 -a 3 )d 1 E 2 + \T 3 31 (a 1 -a 3 )(5E 1 + AE; 



3/ 



1 T 2 21 (E 1 -E 2 )(5a 1 + Aa 2 ) , (45) 

= (E 3 -E 1 )d 2 a 1 + (a 3 -a 1 )d 2 E 3 + iT\ 2 (a 2 -a 1 )(5E 2 + AE 1 ) 

- \V\ 2 {E 2 -E 3 ){ba 2 + Aa 3 ) , (46) 

= (E 1 -E 2 )d 3 a 2 + (a 1 -a 2 )d 3 E 1 + lT 2 23 (a 3 -a 2 )(5E 3 + AE 2 ) 

- |r 1 13 (^3-^i)(5(T3 + 4a 1 ). (47) 

(It is not obvious how the curvature terms in (38) cancel in the transition to 
the tetrad forms; this is shown in Appendix B). 



6 Specializations 
6.1 Type O 

For this class, E v = 0, all constraint equations are trivially satisfied except 
for the Friedmann equation which controls the dynamics, and the evolution 
is that of the FRW models. 
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6.2 Type D 

Without loss of generality we set Ei 
properties are valid: 



E 2 = E ^ E 3 . The following tetrad 



1. From equations (|37| ) and ( |32"| , |33"D we obtain: 

r 1 ™ = r 2 



r 1 



23 



r 2 



23 ) 
13 - 







2. From either fl26|) or ( |2T|) we write: 

(cr 2 - cr^r^a = , 
which has the following two subcases. 



(48) 
(49) 



(50) 



(a) For o"2 = <J\ = o ^ cr 3 the tetrad is free by a rotation in the 
ei,e2 plane. As pointed out earlier, a rotation in that plane can 
be performed so that equation fli~6|) remains valid. This leaves the 
tetrad arbitrary by an initial rotation in a surface x° = const. To 
determine the tetrad completely we perform a further rotation of 
ei , e 2 which preserves (0), where the value of d<p/dx 3 is deter- 
mined from the requirement that r 1 32 = in a surface x° = const. 
From the Jacobi identities (I 



89|) in appendix it follows that 



r 1 



32 







(51) 



everywhere. The tetrad vectors can hence be chosen to be hy- 
persurface orthogonal. This result was also obtained in [TJ. To 
complete the consistency analysis we use the above properties as 
follows: 



i. Constraints 



and ( |35| , p6|) are written respectively as 



d x E -- 



31. 



-ET 3 



31, 



d 2 a - 
d 2 E 



32 



-ET 3 



32 



(52) 
(53) 



with constraints fl3"T| ) and ( P?| ) identically satisfied due to 
equation ([48|). 
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ii. First we recall that the time propagation of constraints (f26 



31Q are identically satisfied. Also the time propagation equa- 
tions (|5S]-(|T]) and are identically satisfied in this class. 
The identity = from fl47|) follows as expected from the 



= in (|37D . On the other hand time propagation equations 



( |45| , ^6|) take the forms; 

= -2Ed 1( r + adiE - aET 3 31 , (54) 
= Ed 2 a + 2ad 2 E + aET 3 32 (55) 

which are identically satisfied on using (|52|) and (^). 
For this class integrability conditions are consistently satisfied, 
(b) For a 2 ^ai. 

We first note that if o 2 ^ o\ = a 3 then from ( |TT| ) we get Ei = E 3 . 
Similarly <J\ ^ u 2 = cr 3 implies E 2 = E 3 . Both cases falls off the 
specified Type D class. So for non- vanishing shear the eigenvalues 
are distinct. Furthermore from the time propagation equation 
(ED we get 

= E 2 -E 1 = (a 2 -a 1 )(3E-^p) (56) 
from which it follows that 

E = \p. (57) 
Taking the time derivative of ([57D gives 

E(a x + a 2 ) = } (58) 
from we get E = (iff a = 0) or 

(<7i + a 2 ) = . (59) 

By a series of three further time derivative of (|)9|) one may show 
that the eigenvalues Ei of the Weyl tensor vanish (i.e., it is type 
O rather than type D), and hence is a FRW solution. 

This proves that for irrotational dust with a purely electric type Weyl tensor 
that is degenerate the shear is also degenerate in the same plane; furthermore, 
the integrability conditions are satisfied. 
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6.3 Type 1 

In this case, E\ E 2 ^ E 3 ^ E\. We deduce the following properties; 

1. From the propagation equations ( |TTD it follows that the shear eigenval- 
ues are also distinct. For if say a x = cr 2 then from ( |TT| ) we get E x = E 2 
which contradicts the requirements of this class. 

2. The tetrad vectors are uniquely determined. 



3. From (|T2| ) we obtain T^o^ = for fi ^ v and hence the spatial tetrad 
vectors are Fermi propagated. 

4. From the new constraint (|39]-fIT]), that is the time development of the 



constraint equations (|32]-|4]), we note that the spatial tetrad vectors 



are hypersurface orthogonal |I| i.e., r 1 2 3 = T 2 3 i = r 3 i 2 = 0. 
For this class further analysis of constraints (|45|-f|7[) is performed below. 

6.3.1 Equivalent Constraints 

We now concentrate on the new constraint (|46|) . If we write the second term 
on the right hand side of constraint (f5|) as 

fa - a l )d 2 E 3 = fa - <Jl)d 2 {E 3 - #i) + fa - o 1 )d 2 E l (60) 

and use (|36|) we obtain the constraint 



(E 3 - £l)9 2 CTi + ((T3 - (7i)5 2 £l 

= |r 3 32 (^ 2 - a 3 )(E 2 - E 3 ) - |r 1 12 (a 2 - a 1 ){E 2 - E x ) 

- Y\ 2 [fa - at)(E 2 - E x ) + fa - (7a) (£ 3 - E x )\ . (61) 

Similar manipulations of the first term on the r.h.s. of (f46|) yields the follow- 
ing constraint: 



(E 3 - E x )d 2 a 3 + (cr 3 - a x )d 2 E 3 
= §r 3 32 (a 2 -a 3 )(£ 2 -£ 3 ) 

- T% 2 [fa-G l ){E 2 -E 3 ) + fa-a 3 ){E 3 -E 1 )]. (62) 



lr 3 32 ((x 2 - a 3 ){E 2 - E 3 ) - lV\ 2 fa - a x ){E 2 - E x ) 

^3 



Each of equations (pl] , p2|) is equivalent to the constraint fl46|), if we assume 



the remaining constraint equations to be satisfied. 
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6.3.2 Consistency test 

To test consistency of the two constraints ( ]6"1"|J6"2| ) we simplify as follows: Let 



X x = d 2 (Ti; X 3 = d 2 a 3 ; Y x = d 2 Er, Y 3 = d 2 E 3 

a=(E 3 -E 1 ); b = (a 3 - a x ) . 
So now ([Jl],^) become 

aX x + bYt = Rt , 

aX 3 + bY 3 = R 3 , 



(63) 
(64) 



(65) 
(66) 



with Ri and R 3 given by the respective r.h.s. Two further constraints are 



obtained from 



by subtracting and adding respectively i.e. 

a(X 3 -X 1 ) + b(Y 3 -Y 1 ) = R 3 -Rx, 
a{X 3 + X l )+b{Y 3 + Y l ) = R 3 + R x . 



(67) 
(68) 

Constraint (|67|j68D are due to the time propagation of the constraint (|36|). 

Next we modify the original constraints as follows: multiplying ( p0|) by 
a = (E 3 - E x ) and © by b = (a 3 - ax) gives: 

(69) 



aX 3 
bY 3 



aXi + Li , 
bYx + L 3 , 



(70) 



where now Li = a x (rhs of (|30"D) and L 2 = b x (rhs of ([36])). Taking a hint 
from QB7| , |65| ) we use constraint ( |69] , [7D| ) to formulate the following two further 
constraints as; 

a(X 3 -X 1 )+b(Y 3 -Y 1 ) = L, + L 3 , (71) 
a(X 3 + X 1 )+b(Y 3 + Y 1 ) = 2aX 1 + 2bY 1 + L 1 + L 3 , (72) 

which have not been time propagated. If we substitute ( ]67|j68| ) into (|69|), ([T0|) 
respectively we obtain: 

L 1 + L 3 = R 3 - R x (73) 



and 



2aX l + 2bY x + L 1 + L 3 = R 3 + R 1 



(74) 
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which is a further set of constraints that are true if and only if (61,62) hold. 
However equation ([F4]) may be reduced to equation fl73|) on using (|65|). This 
reduction is expected since the two equations ( |5"T| , |5"ED result from one con- 
straint, namely ([46]). Thus all we did here was to split one equation (|46|) into 
two equivalent new constraints ( |69l , [70D . And by inserting original constraints 
relevantly modified as fl71), [72|) into the (|67| , |6"8"D and simplifying we obtain one 
new constraint. 



More explicitly the new constraint (^) can be written as 



(E 3 — Ei) T 3 32 (a 2 
+ (0-3 - 0-1 



0-3) 



r 3 32(-^2 



lr 3 32 (a 2 



r 1 i 2 (o"2 - &i) 
e 3 ) — r 1 i 2 (£ , 2 — Ei) 

Ex) - \Y\ 2 {a 2 -ai){E 2 - Ei) 



Vz){E 2 — &3) — 3J- "12^2 - &l){-&2 

- T\ 2 \{a 3 - oi){E 2 - Ei) + (<r 2 - t7i)(S 3 - Ex)] 

- |r 3 32 (ff2 - cx 3 )(E 2 - E 3 ) + |r 1 12 (a 2 - ax){E 2 - E x 
+ r 3 32 [(a 3 - ai)(E 2 - E 3 ) + (a 2 - a 3 )(E 3 - Ex)} 



(75) 



and reduces easily to the identity = 0; which thus shows that (46), the 
basis on which these equations were derived, must be identically true (in 
virtue of all the other constraints that hold). A similar approach for the 
constraints (f45|) and (|47|) yields similar identities = 0. And hence the set 
of integrability conditions (45-47) are also satisfied for Type I fields. 



7 Coordinates 

The vectors of the eigen- tetrad are hypersurface orthogonal for both type D 
and type I spacetimes. So now using the tetrad -coordinate relations (|l4l,|T5|) 
there exists a coordinate system in which line element has diagonal form; 

ds 2 = -V 2 dt 2 + A 2 dx 2 + B 2 dy 2 + C 2 dz 2 (76) 

where A,B,C,V are functions of the spacetime variables (Barnes and Rowl- 
ingson, [[[]). For type D spacetimes, a suitable rescaling of the t coordinate 
to set V — 1 and use of (48) and (52) gives 

ds 2 = -dt 2 + A 2 {x l ){dx 2 + k{x,y) 2 dy 2 ) + a~ 2 {x 3 )g 2 {z,t)dz 2 (77) 
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where a 2 (x l ) is the magnitude of the shear; this has the form assumed by 
Szekeres 0. 



8 Conclusion 

The result that the tetrad vectors are hypersurface orthogonal appears as 
an integrability condition from the H constraint with \i — v, and for the 
given tetrad choice this condition is satisfied. This has been emphasized in 
Ujfl. We have managed to show here that the remaining H constraints for 
fi ^ v are also satisfied without imposing any further geometric requirement 
on the tetrad. We have also checked that the time derivatives of all the other 
constraints are identically satsified in view of the set of constraint equations 
that hold under our given conditions. 

Thus all the constraints resulting from the vanishing of H a b can be con- 
sistently satisfied in solutions with p = = u> a b- More specifically, H a b = 
is equivalent to (6); the time derivative of this equation gives (10); if both of 
these equations are satisfied, there are no further consistency conditions to 
be satisfied resulting from H" ab = or higher time derivatives. 

The latter result was assumed to be valid in the papers by Mataresse et al 
U and Bruni et al , where on neglecting the effect of surrounding matter, 
as carried by H a b (in the absence of a pressure gradient), the model evolved 
as a silent universe. 

We can look at this from the viewpoint of initial data: for initial condi- 
tions that satisfy the constraint equation H a b = (that is, (6)) at some initial 
time t = t , together with the constraint equation H ab = (that is, (10)) at 
that time, then for later times (within the Cauchy development of this initial 
data) the magnetic part of the Weyl tensor will remain zero. Thus vanishing 
of H a b is a consistent condition in an open domain. It is a a dynamical re- 
striction which will remain valid through the evolution of the system; thus for 
example the collapse of fluid elements to a curvature singularity in the form 
of a Zeldovich pancake can be consistently represented by such solutions. 
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A Jacobi Identities 

Following [3], the commutators [e a ,eb] can be written as: 

[e , e b ] = i c abe c , a,6,c = 0...3, (78) 

where the structure constants 7 c a & are related to the Ricci coefficients T c ab 
by: 

l c ab = T c ab - T c ba (79) 

and inversely by: 

Tafcc = \{labc + lcab ~ Ibca) ■ (80) 

The Jacobi identity 

[e b , [e c , e d \\ + [e d , [e b , e c ]] + [e c , [e d , e b \\ = (81) 
can then be written in terms of ^ab as: 

^7^] + 7V7 f % = . (82) 

That is: 

( 123 ) : 5l 7 1 32-'927 1 31 + '937 1 21 = -7 2 327 1 21-7 3 237 1 31+7 1 32(7 2 12 - 7 3 13) , 

(83) 



123 J : ~^ 1 7 2 32+<9 2 7 2 13-537 2 12 = -7 1 3l7 2 12+7 3 137 2 32+7 2 13(7 3 23+7 1 2l) , 

(84) 
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3 
123 



— <9l7 23 + <927 13 — $37 12 = 7 217 13 ~ 7 127 23—7 12(7 31+7 32) , 

(85) 



\ O/JLU 
' 1 N 

v 02 3 / 

2 ) 
031 J 

3 ) 

012 



<9o7 1 32 = 7*32 (#1 -62- 3 ) , 

<9o7 2 i3 = 7 2 i3(-#i + 2 - 9 3 ) , 
<9o7 3 i2 =7 3 i2(-^i -^2 + ^3) • 



(86) 
(87) 
(88) 
(89) 



B Sample Calculations 



We consider here the time propagation of the H constraint. Starting from 
the covariant form [see equation (10)] 

H : hKv j)klm UkE n l;m = (90) 

and using 

{E n l , m ) = (M);m - E n , p u p , m - R n mpq EW + R q mp iE n q u p , (91) 
we obtain as the time derivative of (90), 
= /^)^ . (92) 

The two terms involving the Riemann tensor R a bcd may be shown to be 
zero as follows: 

h(\ V mm Uk [R q mp iE n q u p - R n mpq E q lU p ] 



rf klm u k \R? mpl E\uV - R\ mpq ESu?\ 
+\r, iklm u k [Ri mpl E\v? - Ri mvq E\u?\ 
\ (Term 1 + Term 2) . 



(93) 
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where 



Terml = ^ klm u k [Ri mpl E\vP - R^E^y?] 

Term 2 = \if klm u k \R? mpl E\vP - R mvq E\u?\ . (94) 



Consider now 



Terml = \^ klm u k \R^ mpl E\u^ - R mvq E\uP\ 

= T1 + T2 (95) 



where we set 



Tl = rf klm E\R mpq u k vP ■ (96) 
T2 = ^ klm E\Rl pl u k ^ . (97) 

We write the Riemann tensor R smp i in the terms of the Weyl tensor as 



Rsmpl C sm pl -\- 2 (OspRlm 9slRpm 9mpRsl 9mlRsp) 

R 

--^{g Sp gim - 9sl9pm) (98) 

where 

Csmpl = (VsmabVplcd + 9 smab9plcd)U a U C E , (99) 

and 

g S mab = 9sa9mb ~ 9sb9raa- (100) 

Calculation of Tl: Using (^) in ( p6|) we get 

Tl = ^ klm ESR l mpq u p u k 
= ^ klm ESu p u k [C l 

mpq 2^9 pHqm £ 1 qRpm 9mpRq 9mqRp) 

= Tla + Tlb (101) 



p9qm 9 q9pm) 



where 
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Tla = ■rf 1Am E^ivPu k C i nvq 

= E q l E u u a u c u p u k [rj iklm 7i i 7 ^r} pqcd + ^ klm {g\g mh - g\g ma )g vqcd \ 

= E\E hd ^ klm rf mahnpqcd u a u c u p u k + E\E d m ^ klm uWu k g pqcd 

+E q l E^ klm u m u c u p u k g pqcd 

= E q iE d m r] jklm u i u c u p u k (g pc g q d - 9pd9 q c) 

= E li E d m rf klm u i u k u p u p - W'^Vtti 

= (Symmetry properties) (102) 



and 



Tib = ^ klm E\u p u k [\{g\R qm + g\R pm - g mp R\ + g mq R 



For dust 



rf^Uki^E^R^u* + E\R pm u p ) - \{E\R\u m - E ml R p u p ) 

-?-(E ml u l -E\u m )}. (103) 
o 



Rpm = \pih pm + U p U m ). (104) 



and hence (103) becomes 



Tib = -\ P E\^ klm u m u k -^{E ml u l -E\u m )^ klm u k 

= (Symmetry properties) . (105) 

So now Tl = Tla + Tib = 0. 

Calculation of T2: Using ( |98f) in ( j97D we get 



T2 = jf Um E\R? mvl u k u p 

= ^ klm E\u k u p [C q mpl + \{g%Ri m + g^Rpm) - \{g mp R q i - g m iR q P ) 

-\{g q P 9lm- 9 q i9pm)} 

= T2a + T2b (106) 
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and as in ( |101| ) we set 

T2a = 

= V 



^ klm E\u k u p C q mpl 

■ jklm E i q E bd u a u c u k u p [r] q mab r] plcd + (g q a g mb - g q b g m a)g P icd] 



V jklm V g m ab V P ic d E%E bd u a u c u k u p - rf klm '>E\E ha l u c v? } u m u k (g pc g ld 
(Symmetry properties.) 



9pd,gic) 
(107) 



and 



T2b = r? klm E\u k u v \\{g\R [m + g^Rpm) - \{g mp R q i - g m iR q P/ 
R 

--^{g q P 9im - g q ig P m)] 

= (Symmetry properties) . 
Here also T2 = T2a + T2b = and hence from ( |9~5"D 

Term 1 = 0. 

A similar program of calculations with i and j interchanged 
Term 2 = and the equation ( p2|) takes the final form 

- E n , v p 



= hKr] j)klm u k [(E n l 
In tetrad form ( |11U| ) becomes 



(108) 

(109) 
gives 

(110) 



^(r pe E a - 



dpE T a — 6pdpE T a . p + r r p e E e a — r £ p a E 7 



r e e t 

1 P a- LJ e 



Setting jj, = v in (|111|) yields: 



= rf^uS* 



/3a 



(E a — E^) — 9p(E a — E p 



(111) 



(112) 



If we substitute the propagation equation (0) in (112) and consider the 
values of /i = 1, 2, 3 we obtain the equations 



= ^23(^3-^1)^2-^3), 
= T 2 31 (E 1 -E 2 )(a 1 -a 3 ) , 
= T\ 2 (E 2 -E 3 )(a 1 -a 2 ) . 



(113) 
(114) 
(115) 
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which are the constraints (39-41). 



For values /i^, say // = 1, v = 3, (|111|) becomes 



= d 2 {E 3 — £a) — 2 d2(E3 — E\) + r i 2 (i? 2 — E\) — 9i(E 2 — £a 

- r 3 32 [(E 2 -E 3 ) -6s(E 2 -E 3 )] (116) 

If we substitute the propagation equation ( |i~l~D into (|116|) and use the 
original constraint equations repeatedly we obtain 

= (E 3 - E l )d 2 a l ^ ((7 3 -a 1 )d 2 £ 3 + ir 1 12 (a 2 -a 1 )(5£ 2 + 4£ 1 ) 

- §r 3 32 (£ 2 -£ 3 )(5<7 2 + 4(73); (117) 
The remaining set of fi — 2, v — 3 and // = 1, v = 2 yield 

= (E 2 -E 3 )d l a 3 + (a 2 -a 3 )d 1 E 2 + lT 3 31 (a l -a 3 )(5E 1 + AE 3 ) 

- |r 2 21 (E 1 -E 2 )(5(7 1 + 4a 2 ); (118) 

and 

= (E 1 -E 2 )d 3 a 2 + (a 1 -a 2 )d 3 E 1 + lT 2 23 (a 3 -a 2 )(5E 3 + AE 2 ) 

- |r 1 13 (^3-^l)(5(73 + 4 ( 7 1 ) (119) 

respectively, which are the constraints (45-47). 
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